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ABSTRACT 

Adomian Decomposition Method has been applied to obtain approximate solution to a wide class of ordinary 
and partial differential equation problems arising from Physics, Chemistry, Biology and Engineering. In this 
paper, a numerical solution of delay differential Equations (DDE) based on the Adomian Decomposition 
Method (ADM) is presented. The solutions obtained were without discretization nor linearization. Example 
problems were solved for demonstration. 
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I. INTRODUCTION 

Problems arising from modelling real life problems usually end up with differential, partial differential or 
delay differential equations, etc. Solution of such problems by the usual procedures of mathematical analysis 
such as linearization and perturbation techniques, assumptions of weak nonlinearity etc, necessarily have to 
make slight changes to such problems in order to make them mathematically solvable. Thus, solutions derived 
from such procedures (i.e. linearization, perturbation etc) may not actually reflect the actual behaviour of the 
original problem thereby putting the use of such results to serious limitations. The avoidance of these deviations 
so that physically correct solutions can be obtained is highly desirable. This is the sole objective of the 
Adomiam Decomposition Method (ADM) [1] 

ADM seeks to make possible physically realistic solutions of complex real life problems without the usual 
modeling and mathematical compromises to achieve results. It has been shown [3] that nonlinear systems can be 
extremely sensitive to small changes, thus the ADM is the most desirable option when solving such problems. 
Delay Differential Equations (DDEs) generally provide the best and sometimes the only realistic simulation of 
observed phenomena. [Baker and Paul (1992)] and [Zennaro(1988)]. 

DDE of the first order has the form: y 1 (t) = f (t, y(t), y(g(t)), t > 0, 

y(t) = 0 (t), t < t 0 (1) 
or equivalently as: y (t) = f(t, y(t), y(t — t^), y(t — t 2 ), ... . y(t — r n )) 

where the delays or lags x x are always nonnegative, may be positive constants, 0 < t 1 < r 2 < ••• . T k or 
functions of t , T t = T t (t) or functions of t and y itself, x { = Ti(t,y(t)). Thus essentially DDE types are classified 
according to the complexity of the phenomenon namely 

(i) DDE with constant delay, (ii) DDE with variable or time dependent delay and (iii) DDE with state 
dependent delay. 

Though Delay Differential Equations (DDE ) arise in models through the sciences, DDE with constant delays is 
a large class, especially popular for Biological models [2]. Models in Immunology, Spread of Infection, 
Predator-Prey etc are common instances in which DDEs occur frequently. 

There are lots of differences in finding the solution of ODE initial value problem and a Delay Differential 
equation. For instance given a DDE : 

y 7 (t) = f(t,y(t), y(t- T) ) ; t >t 0 (2) 

y(t) = (p(t), t<t 0 

and ODE IVP: y (t) = f(t, y) , y(a) = c, t > a (3) 
The most obvious difference in solving ODE initial value problem(3) and DDE (2) is that for some t > t 0 , it 
can be that t- T < t 0 ; thus solution of (2) is determined by value at initial point and values prior to the initial 
point, which is called the history . The given initial data must include y(t 0 ), the initial value and values y(t), for t 
< a , say for all t in [a- T, a] . 
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Adding a small delay term in an ODE can change the qualitative behaviour of solutions of the equations. [2], [3]. 
For instance, in the equation 

y^y 2 ©, y(0)=l, t>0 

the solution is not defined for finite points t > t* > 0. But the corresponding DDE 

y 7 (t)= y(t)y(t-l), y(t)=l t<0 has solutions for all t > 0. 
Likewise , all the solutions of the ODE: y 7 (t) = - y(t) are all decaying exponentials, whereas the 

solution of the DDE y 7 (t) = -y(t-7i/2) has solutions of the form y(t) = Asin(t) + Bcos(t). 
For other details of the theories and other qualitative properties of solutions of DDE see [3]. 
The analytic methods for solving Delay Differential Equation include: (i)The Method of Characteristic s,(ii)The 
Bellman's Method of Steps and (iii) The Method of Laplace. 

Until recently the most common numerical method for solving DDE is the continuous extension of Runge 
Kutta methods and its variants [3] . The limitations of these methods have been highlighted earlier in the 
introduction. 

In this paper our focus is on the use of the Adomian Decomposition Method to solve different types of Delay 
Differential Equations. In the next section we present the essential highlights of the method and then go further 
to the implemention. 

II. THE THEORY OF THE METHOD 

Given the deterministic form Fu = g(t) , where F is a nonlinear ODE operator with linear and nonlinear terms 
. This can be written as Lu + Ru + Nu = g 

Here L is the operator for the highest ordered derivative, R is the remainder of the linear term and N is the 
nonlinear term. 

Thus we have Lu = g - Ru-Nu 

Applying the inverse operator L" 1 (the n - fold definite integration operator from 0 to t) 
^>L- 1 Lu = L- 1 g-l7 1 Ru-ir 1 Nu 

d 2 

For example if L — — ^ then L~ x Lu = u — u(0) — tu^(0) 

=^u = u(0) + tu 7 (0) + L~ x g - L~ r Ru - L~ x Nu 

The method show that the solution u(t) can be written as a series thus 

<t) = 2£=o u n =u 0 - L' 1 R 5£ =0 u n - IT 1 X =0 A n (4) 
In the simple case where there is no nonlinear term then the solution u(t) is readily obtained by 
u(t) = X=o u n where the components u n are determined recursively (5) 

In the case where there is nonlinear term Nu , this is defined by the Adomian polynomials 
Nu = YiO A n , where A n are the Adomian Polynomial that can be generated for all forms. 
In [1] there are defined formulas for generating the polynomials. The first four are reproduced below. 

A 0 = /(u 0 ) A t = u t f^(u 0 ) A 2 = u 2 /^(u 0 ) + ^uj fV(u 0 ) 
A 3 = u^Ku^+u^^ul fW(uo)+^i4 / (3) (u 0 ) 

A, = u 4 /^(uo) + [u^^ + u.ui] /(2) (Uo) + U 2 U2 I f^(u 0 )+ut^ f^(u 0 ) 
The Ai s depend only on the components u 0 to Ui 

It has been shown that the series terms approach zero as 
Vmn! ^ m * s tne or( ^ er °f tne highest linear differential operator [1]. 

Another interesting property is that since the series converges rapidly in norm then the partial sum 0 n = 
YJi=o Ui can serve as guide for practical design and lining 0 n = u 

III. EXAMPLE PROBLEMS 

3.1 solve y/(x) = 1 - 2y 2 (j), 0 < x < 1, y (0) = 0, x < 0 

Ly = 1 - 2 y 2 g), => y(x) = L -1 (l - 2y 2 g)) = x - L~W( ~)) 

YoOO = y (0) = x and y n+1 = -2/ Q x A n dx n > 0 (6) 

The nonlinear term is f(y) = y 2 (x/2). The first Adomian polynomial A 0 = f(yo) , then using the generating 
formula in section 2 , we can compute the remaining A n s. 
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A 1 = 2y 0 (|) yi @ f 4 2 = 2y 2 @ y 0 @ + y\ g) 
= 2y 3 (|) y 0 (f) + 2y x (|) y 2 (|) ... etc 



K2 j sv y 2 j sv V2 , s* V2 

Now using (6) we compute the components y n+1 

,7 



^ = J 5^* 6 - 
Jo 



dx = , 

6 

2x 6 x 
dx = — 



yz = fr 



0 24 x ux , 



\2j ~ 3840 



3840 



5040' 



ys 



<D~ 



645120 



Now the partial sum (p 3 is 

1 3 . 1 5 1 7 

yn + yi + y2 + y^ = x — x h — x x 

/i y>3 ^ 19n 5040 



y4 



r 

^0 



48 3840 



2x. 



645120 

- V4 



120 

x 7 dx 



(7) 



362880 



x 3 x^ x 7 x^ 

Now the partial sum (p 4 = X?-i y, = x 1 1 1 

6 12Q 504Q 3 6 2880 



(8) 



Using (8) to compute approximate y(x) values are as shown in table 3.1. These compare favourably with the 
exact solution y(x) = sinx 

Table 3.1 



X 


ADM 


Exact 


Error 


0.2 


0.1986693309 


0.1986693308 


-1E-10 


0.4 


0.3894183422 


0.3894183423 


1E-10 


0.6 


0.5645424735 


0.5646424734 


9.99999E-05 


0.8 


0.717356093 


0.717356090 


-3E-09 











3.2 Solve y / (x) = iefy(0+iy(x); 0 < x < 1, 

yo = i. y«+i = fol el yn(^) + ly n (x)dx, n>o 



y(0) = 1 



(9) 



J"* 1 £ /1\ £ X 

o r 2 + (-)dx = -l+e2 + 



2 

12 3 



1 i 



2 

1 1 



y 2 = -^ + ^4 --e2 + -xe2 --x + -x^ 



4 3 



151 5 1 16 A I 1 I 1 I - 
y 3 = — — + — e2 x + — ei6 x - 6e» x + -xe8 x - -x e2 x + —x 2 e2 x - — — + -e4 x --x 2 + — -x 3 
^ 36 12 9 2 8 32 12x 9 8 48 



5 1, 
— — xe2 

48 



49787 287 i 107 i 512 3 i i 5 i 1 . i 

y 4 = — - + — — e2 x — e* + —-e ^ - 224e32 x + 4x 6?32 x + -xeQ + — x 2 eQ 

^ 216 72 3 27 4 32 
160 3 ii li 151 

+ ^-e 1 6 x - — * e 2 + ^77T X e 2 + ^r* 
72 64 384 72 

This computation continued until y 12 term before obtaining convergence to the exact solution. 

3.3 Solve = + yQ-^ 2 -2 g = ^y (x) + y g) - x 2 - 2, 0<x< 



1, y (o) = 0, y"(0) = 0 ,x < 0 

r x 3 
\ 

rX 3 



yoOO = JT 0 4 ( x ) + y( x ) - x 2 + 2 dxdx evaluated at x=0 gives 

y 0 W = JJ* ;yC0 + yW - * 2 + 2^^//;^y(o) + y (o) - x 2 + 2dxdx 

=> yoW = * 2 - % — , and y n+ i(x) = JJ* ^y n 00 + y n 00 - x 2 + 2dxdxn > 0 



(10) 
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f x f x 3 / . 1 .\ x 2 x 2 

y^ = \ — I x — -—x + — — ——dx dx = 

yx J 0 J 0 4 V 12/4 192 

3/2 = [ I \ (2^4 X4 " 480 X6 ) 

= r r - ( 2483 6 7 

3/3 J 0 J 0 4 V1105920 X 245760 X / + 



191 
2304 

1 



x 4 -- 



36864 30720 
7 „\ 2483 



480 

x 6 dx dx 



70778880 



62914560 



2483 
1105920 X 

x 8 dx dx 



245760 



17381 
566231040 



1351 



5662310400 



AO 



3354533 



74 = 



x 10 -■ 



1038919 



A2 



y$ = 



y 6 



x 12 -■ 



13045963161600 

2579635877 
1763396748627148800 

198433529 
57782984659014411878400 

2438549637881 



765363172147200 



456085441 



A4 



81508116380988211200 
3192598087 



,14 



16 



12820119634666418208768000 X 
12071213366947 



y 7 = 



A8 



227211940956790109811769344000 1977652734087901115801640370176000' 



y (x) = Yit=o Vi gives the partial sum approximation of the solution values. The comparison with exact 
values is given in table 3.2 

Table 3.2 



X 


y(x) ADM 


y(x) exact solution 


0.0 


0.0 


0.0 


0.2 


0.03999993153 


0.04 


0.4 


0.1599895534 


0.16 


0.6 


0.3599513389 


0.36 


0.8 


0.6398650251 


0.64 


1.0 


0.9997300599 


1.0 



3.4 Solve ^ = [1 + y (x)]y (|) , 1 < x < 4, y(x) = 1, 0 < x < 1 
yo = + y00)y (I) dx = 2dx = 2x - 2, 

y ° (I) = x ~ 2 ' yn+i = -Ct 1 + ynW\y (I) dx 

1 2 , 5 9 
yi = 2x - - + -x J - -x z 



359 1 



y2 = -^7 x - 



■ + ■ 



36 2016 126 



5 a 23 

Y b -I 

48 48' 



91 

96 : 



x 7 — — -x 6 + — -x 5 - — x 4 4- -— -x 3 + -x 2 



91 
144' 



1 

4" 



(11) 



ys = - 



594863 4789 , 1565 7 785543 
■x - — — — — x - ■ 



X + ■ 



4064256 870912 



96768 



130056192 



x 7 + 



62221 
6193152' 



, 107329 c 905 4 , 52760122091 

H — — x b — X 4 + ■ 



6635520 



49152 



334764638208 43352064 



393037 o , 891773 q 

X b H —x^ 



83607552 
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883 1 985 47 41287 

x 14 + x 15 + x 13 - x 12 + x 11 



172032 30481920 40255488 196608 29196288 

20445853666660749665209 361507333262608391 g 

3/4 ~ 112067362994533133451264* ~ 326536605461926379520* 



1724824761919 3 130954592996581285 
" 340142297356173312 X ~ 1360569189424693248" 



1415686669370697601 ? 5883282668242351 
108845535153975459840 * + 1332802471273168896" 



13250164149886973 - 1744394411483495 , 
x + x 



2591560360808939520 58101081182011392 
22977368275584935 n 1605814479948579829 



x a + — - x 1Q 



15832077840578248704 979609816385779138560 



618262972935793 „ 3984471945937759 
x 14 + x 15 



10719636037891522560 148425729755421081600 



2740702697436551 _ 84668516381708137 



-x 13 + x 12 



173264321265511956480 293216543680097157120 ' 



189594690223731469 „ 16517253935503343176150636474441 



213803729766737510400 187115261819515399230021697536000 



2396541839181487 _ 1 

ylO _| ' 

373184691956487290880 28803570853478400 ' 



1 on 48389 



-x 30 + x 29 



330363536670720 398550570639556608 " 



90983 nn 1628218573 



-x 28 + x 27 



30540273612226560 32487060237283491840 

18900269 9A 26357624082151 

x 26 + x 25 



30809588194344960 4664808649456091136000 ' 
11578815968243 n4 1995670980206023 



-x 24 + x 23 



290254760410601226240 9179306797985263779840 ' 

121196867155951 99 4705550936141747 
x 22 + - -— — x 21 



133033431854858895360 1632683027309631897600 ' 



103272035359 9n 337671295473833 
x - x 1 ^ 



15670810949713920 34531700392077557760 ' 



144654049387121 1Q 10333140079326283 

x 18 - x 17 



2332404324728055680 1850374097617582817280 
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With the aid of Maple software y 5 , y 6 , y 7 , were calculated . The approximator y(x) = 2^=q ¥i gives 
approximated solutions which converges to the exact solution. 

III. CONCLUSION 

In this paper, a numerical method based on the Adomian decomposition method (ADM) was 
implemented to solve both linear and nonlinear delay differential equations. The numerical results obtained by 
using the Adomian decomposition method compares favourably with the exact solutions. The attractive 
property of the method is that it is implemented directly in a straightforward manner without using restrictive 
assumptions, linearization, and discretization. The efficiency of the decomposition method gives it much wider 
applicability which has to be explored further. 
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